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A phenomenological theory for the behavior of 
axial rotation of main-sequence stars is proposed he 

jv 

by considering the effect of braking. It does not 
specify the physical process of braking but does . 


provide a statistical model by Vv^ich histograms of ' 
observed rotational velocities of stars of different 
spectral types can be explained in terms of the braking 
strength. It shows that rotation of all stars including 
those of 0 and B type has been braked in various degrees 
during their course of evolution. 

As a result of braking, the stellar angular 
momentum is transferred outwards to the surrounding 
nebula that may be regarded as the remnant of star 
formation. This leads plausibly . to the formation around 
the star of a planetary system from the nebula. If the 
angular momentum is not further dissipated from the 
nebula, we can estimate the size of the planetary system 
and its probability of occurrence. It shows that the 
size of a planetary system around any star is critically 
dependent upon the mass in the nebula as compared to 
that of the star itself. 
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Axial rotation of the rnin-sco 
abruptly at about F5 (Struts 1030). 
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attributed to the braking effect of various mechanisms reviewed 
in a previous paper (Huang and Struve 1960). A more efficient 
braking mechanism than those previously proposed has recently 
been suggested by Schatzman (1962) as due to mass loss through 
stellar magnetic activities. 

In the present paper we do not intend to discuss the 
physical process of how the braking of stellar rotation is 
actually Effected. Rather, we will show phenomenologically 
that the observed distribution of rotational velocities of 
stars belonging to different spectral types can indeed be 
understood in terms of braking strengths. 

As the result of a statistical study (Struve 1945; Huang 
and Struve 1954) it has been concluded that the rotational 
axes show no preference in their directions in space. This 
indicates that the rotating stars have not acquired their 
angular momentum from 


the galactic rotation of the pre- 
stellar gaseous medium. Instead, their angular momenta must 
have been derived from a random process, such as collisions 
of interstellar clouds, turbulence, etc. (Huang and Struve 
1954). Consequently, we may expect that the angular momenta 
per unit mass, h, of stars are originally distributed according 
to the Maxwellian law. Finally we would like to call attention 
here to the fact that this observational consequence has not 
been fully appreciated. Many theoreticians still take the 
rotational angular momenta of stars as due ultimately to the 
galactic rotation of the tenuous pre-stellar medium (e.g. 
Edgeworth 1S46; Hoyle I960). 

For a given spectral type the geometrical radius and the 






radius of gyration do not vary 
the equatorial rotational velocity, v, to bo 
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proportion to h. If we set J/C — ly '/u ]K > where lT )n denotes the 
most probable equatorial rotational velocity, the original 
distribution of x in a given spectral range should be given by 





according to our assumption. Here Z 0 is the upper limit of 
rotational velocities because of rotational instability and 
is a numerical factor to normalize the function j 
if we assume that the star rotates as a rigid body, Z is 
obviously related to the angular momentum per unit mass, h, 
of the star as follows: 


fZj . Now 





where R is the geometrical radius and 1-Jc the radius of gyration 

j} _p 

of the star. We may define ;;• = by equation (2) when X - ^ . 

r, p 

Then interstellar clouds with f ) 7 ft could very possibly lead to 
the formation of close binary systems. Since at present we 
are interested mainly in the behavior of the observe- 5 distribution' 
at its small- velocity end which is not significant ..fleeted 
by the high-velocity tail, we shall set, for simplicity of 
calculation, and consequently 'Q — ( . 

As a simple model, we propose that the effect of braking 
is to reduce h of all stars in a given spectral range by a 
constant amount, h, . Define j = X, when is substituted 

in equation (2). It is then obvious that Zj measures the degree 
of braking of axial rotation of the stars within a spectral 
type, namely, the stronger the braking mechanism, the greater 
the value of Xj . In this way we are able to obtain the 
distribution of rotational velocities after the braking mechanism 
has been applied. 

Needless to say, this statistical model is very idealized 

Z 

and represents an extreme case. Most likely is not a constant 
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.even within a spectral sab-typo but increases with the original 
value of h. However, since the observed distribution is not 
well determined, we do not see any advantage at this moment 
to build an elaborate model which us* ar_y would involve many • 
adjustable parameters. On the other ..and, because of this 
oversimplified model, the quantitative results obtained in 
this paper represent only a very rough estimate. - 

If we accept this simplified model, the distribution of 
equatorial rotational velocities after braking has taken its 
toll becomes V' 
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and ^(y- 0 ) denotes the delta function. 

If denotes the distribution of observed velocities 

y, we have (Kuiper 1935; Chandrasekhar and Munch 1950) 

r4 . ; 


>*i)= 'd 


j 

'll 

Substituting equation (3) into equation (5) we ooxain 
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The observed distributions are given in the form of 
histograms, i.e., 
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In Figure 1 we have shown 11 computed histograms of 

§ (J, *t) according to equations (5)-(7) for // values of X { 

(0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 1.0, 1.4, 1.8, 2,6) with<d^-0*2. 
This sequence of histograms illustrates the effect of the ^ 
braking strength on the distribution of observed rotational 
velocities of stars in a spectral type. The case ~ O 

corresponds to that of no braking, while the case = 2.6 leads 

practically to a complete stop of axial rotation by braking. 

These diagrams should be compared with the observed 
histograms. Unfortunately, as we can see in our previous paper 
(Huang 1953) , the two sets of observed diagrams derived from 
two uncombinable sees of observational data show slight 
differences. The two uncombinable sets of observational data 
are respectively from the three-prism and two-prism spectrograms • 
that have been taken during the years by astronomers at the 
Lick Observatory. There are 1103 stars for which rotational 
velocities were measured from three-prism spectrograms and 
445 stars for which rotational velocities were measured from 
two-prism spectrograms. Therefore, the histograms of rotational 
velocities derived from the three- prism spectrograms are 
definitely more reliable than those from the two-prism spectro- 
grams from a consideration of both the dispersion and the number 
of stars studied. For this reason we shall compare the computed 
histograms in Figure 1 with the observed ones obtained from 
three-prism histograms, the latter being reproduced here 
in Figure 2. Finally, it may be noted that there are available 

better observational data for rotational velocities (Slettebak 


We have mentioned before that by assuming ~X -v' >/o we 
introduce an approximation. While this approximation does 
not seriously modify the distribution function at the lower 
end of observed rotational velocities, it makes an appreciable 
difference at the upper end. Thus, when we compare histograms 
(those of r= 0.3, 0.4 and 0.8) in Figure 1 with those in 
Figure 2, we can immediately see a systematic difference 
arising from this approximation, since in all three histograms 
obtained from observed data, the distribution function drops 
much more rapidly than what have been calculated on the 
assumption of jC, — 5> fo . Thus, a cut-off of rotational velocities 

i 

at their upper end is evident in all three cases in Figure 2. 

This shows that rotational instability has played a role not 
only in 0- and 3-type stars but also in F- and perhaps even 
later-type stars sometime during the course of their evolution. 
This result may be significant to our search of clues for 
star formation. 

The strength of braking varies with, the moment of braking 
force and the duration in which braking acts. If indeed the 
braking is due to magnetic activities which are in turn 
induced by convection and rotation itself in the early phase 
of the star’s evolution as Schatzman (1962) has suggested, 
we would expect from the present study that magnetic activities 
have occurred in stars of all masses in their early phases of 
evolution, differing only in the strength and duration of such 
activities. This conclusion is in agreement with Pvoeda T s 
(1964) interpretation of flare stars based on Hayashi’s (1961) 
results of evolution on one hand and Fowler, Greenstein and 
Hoyle’s (1962) theory of formation of light elements — Li, Be, B - 
in the early phase of the solar system on the other. 

Thus far, we have followed the reasoning of both Schatzman 
and Poveda that the magnetic field is generated in the star as 
a result of convective motion and differential rotation. 

However, it may be noted that it would not modify our 
conclusion if we assume that the strong magnetic field is 
derived from accretion of interstellar magnetic lines of forc e 
during the period of star formation. If so, the strong magnetic 
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here. However, our crude measurements of stellar rotational 
velocities provide the only uniform set of data which cover such 
a large number of stars that make the histograms in Figure 2 
statistically significant. 

In Table 1 we give the results of our determination of I 
Xj and for each of three observed histograms after a 
comparison with those in Figure 1, being determined by 
the shape and 1 T ;| by the scales of abscissas of both observed = 
and computed diagrams. It should once more be emphasized that 
these values are very tentative as both theory and observational 
data (especially the conversion of measured widths of the 
spectral line to rotational velocities) are crude. However, 
there is no doubt at all that is least (but -'definitely non- 
vanishing) for 0 and B- type stars, increases only slightly 
for A-type stars and becomes much larger for F0-F5 stars. Since 
few single stars of spectral types later than F5 are found to 
be rotating with a measurable velocity, we must conclude 
that for these late- type stars is greater than 2.6. 

Yttiile we cannot determine x -accurately because of the 
uncertainty in observational data, the general behavior does 
indicate that our simple statistical model for the braking 
effect represents very well the gradual change of the distribution 
of observed rotational velocities from 0 and B through A and 
F0-F5, and finally to later-type stars. Especially it explains 
the long puzzled fact that a large number of 0, 3, A stars 
do not rotate appreciably, although their average rotational 
velocities as group are high. 

The previous result induces us to conclude that contrary 
to our previous belief, braking is not limited only to stars 
of spectral types later than F5 but extends to 0, B, A and 
early F stars as well. It suggests that along the main 
sequence all stars show evidence of having been braked. However, 
the braking mechanism is weak at the upper branch of the main 
sequence and increases its strength rapidly after F0. It 
becomes so strong after F5 that rotation ox all- stars is 
practically stopped by it. 
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activities in the early phase would represent the 

because of convection,, of magnetic energy that is 
there . 


dissipation 
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It has often been suggested (e.g. Huang 1953) that the 
disappearance of rotation after F5 may indicate the emergence- 
of planetary systems because it is difficult to understand how 
the original distribution of angular momenta of the late-type 
stars should be so radically different from that of early ones. 

This is especially true if we follow the reasoning in the 
previous section where the original distribution of angular 
momenta per unit mass of stars is assumed to follow the 
Maxwellian distribution (which may be truncated at the high- 
value end because of rotational instability). Indeed, if 
angular momenta that have been dissipated 

from the stars remain in their neighborhood, formation of 
planetary systems would appear inevitable. Following this line 
of reasoning we shall make in this section an estimate, though 
a very crude one, not only of the probability among single 
stars that will possess planetary systems but also of their 
general behavior with respect to the spectral type of their 
parent stars. 

Before we shall understand why our planetary system was 
formed in a state as it is, we are not expected to derive, 
simply from. a consideration of the angular momentum, the 
detailed structure of planetary systems around stars. In order 
to obtain a general behavior of any planetary system around a 
star without going into the detail structure we may define the 
concept of an "equivalent planet”, which is a fictitious 
planet that would be moving in a circular orbit with an 
angular momentum pex* unit mass, h, equal to the average value 
of the entire mass outside of, but belonging to, the star. 

Thus, in the case of our solar system, we have (Allen 1955) the 
total angular momentum of the planetary system equal to 3.15 x 10^ 
gm cm^/sec , while its total mass is 2.68 x 10^® gm. 


8 


J> f 

Consequently, its h } , denoted hereafter by ^ 
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Then it is easy to find that the equivalent planet of our 
own planetary system is located at a distance 

a s - j. 0 cj-2 x 1 0 1 ^ 




from the sun (i.e. between Jupiter and Saturn as would be 
expected) and revolves with a linear velocity 
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around the sun. It is obvious from our definition of the 
equivalent planet that it does give us a measure of the extent 
of any planetary system although it conveys no idea about its 
total mass and mass distribution in the system. 

Let us now assume that the gaseous remnant around a star 
after its formation is m and the mass, radius and radius of 
gyration of the star itself are respectively M, it, and Rk. 
Therefore, the total angular momentum of the star is 
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After braking, all of its angular momentum is transported 

into the surrounding medium (of mass m) according' to our 

theory if % ^ 7/, X, l^.yand a constant part - (equal to _Cl, ~ Mi'b 

of it is transported out if > l/j . However, we shouad 

realize that the angular momentum in the nebula that will 

evolve to become a planetary system could be different from 

what has been fed^into it by the star, because some angular 

momentum, say f> , may be originally associated wa. an it and 
~ w 0 ^ 
has never gone into the scar. Since both . the scar and 

its surrounding nebula are supposed to have been formed from 

the same interstellar medium, we may reasonably assume that 
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points to the same direction as the stellar 


the vector , ^ L , 
angular momentum does. Hence, we can consider their magnitudes . 
only and the angular momentum associated with the nebu_a (of 

mass m) after the completion of braking process has the following 

* 

value: 
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For stars with slow rotation which indicates smallness of h 
of the pre-stellar medium, must be small and may be neglected. 
However, for stars of rapid rotation and especially for those 
which have passed through the stage of rotational instability 
in the course of their evolution, h of the prestellar medium 
must be large and_fl 0 large too. 

If we include the possibility that a part of the angular 
momentum in the nebula may be dissipated away (for example 
through the loss of mass), JT^may be even negative. Thus, _ L 0 
cannot be estimated from the present consideration. However, 
as we shall see, this uncertainty can be formally circumvented 
in the following calculation. 

Since 
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where a is the radius of the equivalent planet*s orbit around 
a star and 
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for the planetary system of the sun, we may introduce a new 


variable 
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if 0 is positive and zero to > if >'? is negative. 
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Therefore, the lower and upper limit of the radius of the 
equivalent planet *s orbit denoted respectively by aQ and a^ 
are given by 
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It is easy to derive from equations (I) and (13) that the 
distribution functional) with respect to 



If we further let 
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Since from j to ^ it involves only a horizontal translation, 

the shape of the distribution function does not depend upon_/7. • 

jTuo iS positive, we simply translate the entire curve to 

the rignt, thereby increasing systematically the size of the 

expected orbits. If__Q_ is negative, we translate the curve. 

to the left and set the probability to zero everyv/here that 

is on the left side of the origin. In this way we are able 

formally to take into account any value of Ar^. 0 . Physically, 

however, we should remember that O likely varies from case 

r~* < -o 

to case and moreover we do non know how it does vary (say, 
with respect to M, h, m, h* , etc.). Hence, it is one of the 
uncertain factors that prevent us to derive the probability 
of occurrence of a planetary system of a given size. 

Since we can obtain easily the distribution with respect 
to ^ from that with respect to ^ , we shall consider hereafter 
only the latter which, as we see from equation (24), depends 
upon two parameters, y\- and ~f, , or equivalently /\ and ^ because 
of equation (23). While ^ has been determined from observed 
data, 7q involves an unknown parameter, namely m/M. 

From the values of Aj and in Table 1 we can compute 
'_~r \ ) for each of the three cases given there. We shall 
.use the mean radius R/R^ = 3.9 and mean mass M/M^ — 5.9 for 
0 and B stars, R/R^ = 1.74 and - 2 for A stars and R/R^ = 
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and H/Hq = 1.4 for F0-F5 stars. In all cases we adapt 


- 0.05. 
and (23) that 


follows from Table. I and equations (20) 
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lat-./M is equal to C.S x 1C~' for 0 and B start, equal to . 
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Consequently, we may roughly set 
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as representing all stars from 0 to F5. 
(24) , when combined with equations (17) , 
gives 



If = 0 , equation 
(18) and. (23), 



which states that the maximum size of the planetary system 
is directly proportional to the mass of its parent star 
and inversely proportional to the square of its total mass. 

It should be noted that the numerical coefficient in equation 
(24) does not apply to stars later than F5 for which we do 
not have any observed data for making the calculation. 

We have seen how critically does the size of a planetary 
system depend upon the mass in the nebula. A further complication 
is the dissipation of mass in the nebula and accompanying 
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dissipation of angular momentum before the formation of a 
planetary system of finite bodies. The difficulty of this 
point can be seen from the fact that even in our own planetary 
system there is no consensus concerning the amount of mass in 
uhe solar nebula in the beginning, as some investigators 
advocate a value of rr./l.l = 0.1 while others use considerably 
smaller values. 

Since the actual size of any planetary system must be - 
equal to or less than the upper limit, we have for our own 
planetary system f ^ R c . It then follows from equations 
(r.y and (IS) that 



The ratio m/M for the solar system at present is equal to 
—3 

1.3 x 10 . In the early days it must be much greater than this 

value. If we now take iT ; , to be the same as that determined 
in Table 1 (200 Km/sec), we have 
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5.9 x 10‘ 
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Therefore, it follows from equation (2 p) that if _/2 — X > 2.2. 
Equation (1) shows that the probability of having^ >2.2 is rather 
small. This would put our planetary system as an unusual 
case, namely its chance of occurrence is- small. The other alter- 
native is that is appreciable. This . means that a significant 
portion of the angular momentum that is now found in our 
planetary system has never belong to the sun. The present 
consideration cannot decide which alternative is the more 
plausible one. 

Thus, we have seen many uncertainties concerning planetary 
systems around the stars in general.. They cannot be resolved 
for the time being. What we can predict from the present 
simple theory is that the size of a planetary system increases 
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with a decrease of its total mass and that the distribution 
of sizes (measured by y of the equivalent planet *s orbit) 
behaves like a Maxwellian distribution of velocity magnitudes 
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Therefore unless_X2 5 has a large variation among individual 
systems, their sizes do not vary greatly from one to the other. 
Most of them are crowded at the upper limit. The truncation 
of systems around F-type stars occurs at relatively larger 
values. Consequently, the spread of sizes is also larger, 
although there is still an appreciable number of them at the 
upper limit of their size. On the other hand, the sizes of 
planetary systems around late-type (later than F) stars are 
expected to be distributed like a complete Maxwellian curve 
with a wide spread of values. They should not show the 
tendency to accumulate at any particular size like systems 
around early-type stars do. All these consequences would be 
modified if SL Q varies greatly from one system to another; in 
such cases the distribution of planetary sizes can be derived 
only when we know the variation of_ L o* 

Finally, we would like to ask if it is possible that 
there is no mass left after the actual processes of star 
formation (i.e., m = 0). In that case no planet would be 
formed. Indeed from the point of view of star formation we 
do not have any compelling reason to assure us that there must 
be some mass left behind to form planets after the star 
itself is formed. On the other hand, the rotational behavior 
of stars does suggest that there must be some mass left 
behind since axial rotation could not be effectively braked 
without the presence of such mass. around the star, although 
we cannot at present estimate the amount of mass that is there. 
It is therefore the main point of this paper, if nothing else, 
to show the behavior of axial rotation of stars indicates the 
high probability of occurrence of planetary systems around sta: 


The present paper perhaps represents the first serious 
attempt to link the problem of occurrence of planetary syste: 
in the stellar universe with some observat ional facts. The 
uncertainties in both theory and observation force us to aca 
an over-simplified model. Thus, it should be regarded as a 
pioneer exploration of a nebulous field of .Learning rather 
than a legitimate treatment of something that is clearly 
understood. As a result we emphasize only the qualitative 
conclusions that have been derived here. 


It is my pleasure to express my thanks to Mr. Clarence 
T7ade, Jr. for the numerical integrals in equations (7) and 
(8) on the IBM 7094 digital computer at Goddard. 
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Fi gure 1 


Figure 2 


LSGBND3 OF FIGUTiFS 


Computed hi programs of projected rotational- 
velocities of "roups of stars that have suffered - 
different degrees of braking. 7/e have plotted here 
11 histograms for 11 values of j which, according 
to our model, measures the strength of braking of * 
stellar axial rotation. 


- Observed histograms of projected rotational 

velocities of three groups of stars, 0 and B, A, 
and F0-F5, taken from a previous paper. The 
histogram of projected rotational velocities for 
later type stars will be consisted of only a single 
column of small velocities. These histograms 
should be compared with those in Figure 1. The 
trend of increasing braking strengths from 0 and 
B, through A and F0-F5 to later-type stars is 
unmistakable. From the comparison we may make a 
rough estimate of ~jQ for each of these three groups 
of stars. 



